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ABSTRACT 

The astronomer E.V. Pitjeva, by analyzing with the EPM2008 ephemerides 
a large number of planetary observations including also two years (2004-2006) 
of normal points from the Cassini spacecraft, phenomeno logically estimated a 
statistically significant non-zero correction to the usual Newtonian/Einsteinian 
secular precession of the longitude of the perihelion of Saturn, i.e. Atusat = 
-0.006 ± 0.002 " cy-^; the formal, statistical error is 0.0007 " cy-^ It can be 
explained neither by any of the standard classical and general relativistic dy- 
namical effects mismodelled/unmodelled in the force models of the EPM2008 
ephemerides nor by several exotic modifications of gravity recently put forth to 
accommodate certain cosmological/astrophysical observations without resorting 
to dark energy/dark matter. Both independent analyses by other teams of as- 
tronomers and further processing of larger data sets from Cassini will be helpful 
in clarifying the nature and the true existence of the anomalous precession of the 
perihelion of Saturn. 

Subject headings: gravitation — relativity — celestial mechanics — ephemerides — planets 
and satellites: individual (Saturn) 
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Introduction 



At present, the best theory of the gravitational interaction available to us is the 
Einsteinian General Theory of Relativity (GTR) which has passed so far many observational 
tests concerning orbital motions an d propagation of ele c tromagnetic wa ves in the (inner) 



Solar System with excellent results (INi 



2005 



Will 



2006 



Turyshev 



20081 ). Deviations from 



the expected behavior have been detected in the hyperbolic motions of the Pioneer 10/11 
spacecrafts after they passed the threshold of approximately 20 Astronomical Units (AU) 



( lAnderson et al. 



19981 ). but it is unlikely that the Pioneer anomaly ma y be ascribed to 



long-range mo difications of the know n laws of gravitation (jlorio 



fiyby anomaly (lAnderson et al. 



2007ah . The so-called 



20081 ) consists of a small but unexplained increase in the 
velocity of several interplanetary probes (Galileo, NEAR, Rosetta) moving along their 
hyperbolic orbits at their closest approaches to the Earth; at present, no conventional 
explanations in terms of known physics have been found. Another anomalous effect 
which has recently attracted at t ention is the se cular increase of the Astronomical Unit 



( Krasinskv and Brumber 



( iLammerzahl et al. 



20081). 



2004 



Standish 



20041 ). For an overview of such topics see 



In this paper I will focus on a recently detected non-standard feature of the motion of 
Saturn which, if confirmed as a genuine dynamical effect by further, independent analyses, 
may be added to the list of the Solar System anomalies not explained by known mundane 



causes. 



The astronomer E.V. Pitjeva has recently processed a huge data set (1913-2007) of 
planetary observations of various kinds including also three-dimensional normal point 
observations of the Cassini spacecraft (200 4-2006) with the refined dynamical models of 



the EPM2008 ephemerides flPitieva 



2008a|). They encompass also the action of Eris, the 



other 20 largest trans-Neptunian objects (TNOs) and a TNOs' massive ring in addition 
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to the usual Newtonian (solar quadrupole mass moment J2, N-body interactions with the 
major planets, 301 biggest asteroids, massive ring of the small asteroids) and the general 
relativistic Schwarzschild-like forces of order (9(c~^). As a result, she produced a global 
solution in which she phenomenologically estimated, among many other parameters, a 
correction Atusat to the standard Newtonian/Einsteinian secular, i.e. averaged over one 
orbital revolution, precession of the longitude o f the perihelion zu of Saturn whose orbital 



parameters are in Table [TJ It is (jPitjeva 



2008a) 



Awsat = -0.006 ± 0.002 arcseconds century"^ (" cy"^) 



it is not compatible with zero at 3 — a level. Concerning the quoted uncertainty of 0.002 
" cy~^, it is important to note that it i s not the forma l error which is, instead, three times 



smaller and amounts to 0.0007 " cy ^ (jPitjeva 



2008bl ). There should be considere d also 



the po ssibility that the realistic uncertainty may be up to 10 times the formal one (jPitjeva 



2008bl ). but I believe that until no other independent determinations to be compared 
with the one of eq. ([T]) will be available, this cannot be decided. Thus, throughout the 
paper I will rest upon the result of the fit of eq. ([1]). However, in Section [2] 1 will present 
some considerations, based on the action of the trans-Neptunian objects, disfavoring the 
possibility that the real uncertainty can be as large as 0.007 " cy~^. Arijsat takes into 
account, by construction, any unmodelled/mismodelled dynamical effects affecting the orbit 
of Saturn. Previous estim ates based on the EPM2006 ephemerides, which did not include 



the Cassini data, yielded (jPitjeva 



20061 ) 



Azjjsat = -0.92 ± 0.29 " cy^^ (formal error . 



(2) 



Waiting for independent confirmations of eq. ([T]) by other teams of astronomers and 
further data analysis including, hopefully, more Cassini normal points, in this paper I 
will address the following questions, a) May some known standard physical effects, not 
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Table 1: Orbital parameters of Saturn at the epoch J2000 with respect to the mean ecliptic 
and equinox of J2000. a is the semimajor axis in AU, e is the eccentricity, i is the inclination 
of its orbit with respect to the Sun's equator in deg, Pb is its orbital period in yr. 

a (AU) e i (deg) Pb (yr) 

9.53707032 0.05415060 4.6655 29.453 

properly modelled, or unmodelled at all, be the cause of the estimated anomalous retrograde 
precession of Saturn? h) Could some of the recently proposed modified models of gravity, 
not modeled in the EPM2008 ephemerides, account for Awgat? 



2. Possible explanations of the anomalous perihelion precession of Saturn 



In Table [2] 1 quote the analytical expressions and the nominal values of the secular 
perihelion precessions of Saturn due to the known dynamical effects of classical and 
relativistic origin, along with some exo tic forces recently proposed to explain, e.g., the 



dark energy (iDvali et al. 



Pioneer ano malv fjAnderson et al. 



19981 ). the cosmo logical expansion wit hout resorting to 



20001 ) and the fiat rotation curves of galaxies (iMilgrom 



19831) 



without invo king dark matter; also the acti ons of a spherically symmetric di stribution o f 



dark matter (IKhriplovich and Pitjeva 



20061 ) and the cosmological constant flKerr et al. 
20031 ) in the Solar System are considered. It turns out that the majority of the considered 
effects, modelled or not in the EPM2008 routines, cannot explain both the sign and the 
magnitude of Arogat because they induce prograde perihelion precessions. In particular, 
since the modelling of the trans-Neptunian objects (TNOs) is certainly not yet complete, 
it may be argued that the uncertainty in their total mass can induce a mismodelled 
perihelion precession of Saturn large enough to explain eq. ([1]). Although it is reasonable 
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to assume the modeling of the action of the TNOs as still preliminary, it seems difficult 
to attribute the determined anomalous apsidal precession o f Saturn to them because the 



2007bl ). Note that if one assumes 



perihelion precession due to them is prograde (llorid 
for the uncertainty in Awsat the value 0.007 " cy~^, i.e. 10 times the formal, statistical 
error, a prograde extra-precession of +0.001 " cy~^ would be allowed. If attributed to a 
still imperfectly modeled action of the TNOs, such a positive apsidal precession would 
imply 5mTNOs = O.OIM^; it is just the nominal value of the Classical Kuiper Belt Objects 



(C KBOs) which constitut es about 70% of the entire popu 



by 



Bernstein et al. 



n 



ation of the TNOs obtained 



( 20041 ) with non-dynamical techniqueqj. This would be in contrast 
to the fact that Pitjeva did actually model the action of the 20 largest TNOs and of a 
massive ring, i.e. it is difficult to believe that a mismodelling as large as 100% may have 
occurred given the exquisite level of modeling of the EPM ephemerides. The general 
relativistic gravitoelectric Scwharzschild-like force of order (9(c~^) has been modelled in 
EPM2008. Since the resulting perihelion precession ca n be writte n in terms of the standard 



Eddington-Robertson-Schiff PPN parameters (3 and 7 (jWill 

nGM 



19931) as 



with 



2 + 27-/3 



(3) 



(4) 



it may be argued that, in principle, Awsat niay be explained in terms of a deviation Az/ of 
V from its general relativistic value vqtk = 1- However, this would imply 



Ai/ 

which is contradicted by several other 
the Solar System (}NJ 



-0.428 ± 0.142 



(5) 



2005 



Wii: 



2006 



independent de terminations of (3 and 7 throughout 



Turyshev 



20081). 



^They used the ACS camera of the Hubble Space Telescope. 
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Concerning the forces able to induce a negative perihelion precession, the Newtonian 
N-body interactions with the major planets yield the largest retrograde effect; since it 
is mainly due to the Jupiter, the uncertainty in its mass might, in principle, induce a 
mismodelled precession able to accommodate eq. ([1]). In fact, the answer is nega tive 



becau se the mass of Jupiter is presently known with a relative accuracy of 1 x 10"^ (iPitjeva 
2008d ) which yields a mismodelled precession two orders of magnitude smaller than eq. 
([1]). The retrograde perihelion precession of order (9(c~^) due to the general relativistic 



gravi tomagnetic Lense-Thirring force generated by the Sun's rotation (ILense and Thirring 



19181 ). not modeled in the EPM2008 ephemerides, is smaller than A'cusa.t by four orders 



of magnitude. With regards to the putative exotic forces con sidered here, the Dvali- 



Gabadadze-Porrati ( PGP) modified grayitv (jDvah et al.ll2000l ) predicts a negative secular 



perihelion precession (ILue and Starkman 



20031) of -0.0005 " cy-i; it is too small to explain 



eq. ([T]) which rules out its existence at 2.75 — a level. Incidentally, let us note that the 
positive value of the DGP effect, related to the self-accelerated branch of the cosmic 
expansion, is ruled out at 3.25 — a level. The anomalous acceleration experienced by the 



Pioneer 10/11 probes at the Saturn orbit flNieto and Anderson 



20051 ). if attributed to a 



constant and uniform extra-force directed towards the Sun acting also on the planets, 
induces a retrograde secular precession of the perihelion of Saturn four orders of magnitu de 
large r than eq. ([I]). The force quadratic in the radial velocity considered in (jStandish 
20081 ) as a possible explanation of the Pioneer anomaly induces a retrograde precession 



of the perihelion which, for Saturn, is o ne order o: 
incompatible with it at 11 — cr level; in 



determinations of Awsat- Also MOND (IMilgrom 



lorio 



magnitude larger than eq. ([T]) being 



20091) it was still compatible with less recent 



19831 ) causes retrograde secular perihelion 



precessions; however, in the case of Saturn they are either too small or too large with 
respect to eq. by several orders of magnitude. 



It may be argued that some mutual cancelations among different unmod- 
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elled/mismodelled effects may have conspired to yield just the estimated value of 
Aujsat, but an inspection of Table [2] shows that this seems to be a very unlikely possibility. 

The detected anomalous perihelion precession of Saturn may be used to phenomeno- 
logically constrain the existence of an unknown constant and uniform acceleration directed 
towards the Sun continuously existing in the spatial regions swept by the Saturn's orbital 
motion during the time interval spanned by the data set used (1913-2007) covering about 
four orbital revolutions of Saturn. It is 



^sat = -(9.14 ± 3.04) X 10"^^ m s" 



(6) 



However, its existence in the inner regions of the Solar System is ruled out by the estimated 
corrections to the Newtonian- Einsteinian perihelion precessions of Venus, Earth and Mars, 
as shown by Table [31 Concerning eq. (Q, it must be noted that it could not be reproduced 
by a Yukawa-like term 

Ay.k = -^: ^(l + y)e xp (-- ) (7) 



evaluated at the Saturn's orbit. Indeed, in (lorio 



20081 ) it has been shown that the estimated 



corrections to the standard precessions of the perihelia of the inner planets constrain a and 
A to a < 4 X 10"", A < 0.18 AU; such values in eq. ([7]) yield for Saturn 



A 



Yuk 



-2 X 10"^^ m s"^ 



On the other hand, typical values for a and A able t o fit astrophysica 
distant galaxies, i.e. a = -3 x 10"^ A = 33, 000 AU ^Moffat and Toth 
for Saturn 



observations of 



Avuk = 1-98 X 10"^^ m s 



20071), would yield 
(9) 
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3. Discussion and conclusions 

Based on the analysis presented, summarized by Table [2l I conclude that the recently 
estimated anomalous retrograde apsidal precession of Saturn Awsat = —0.006 ±0.002 " cy~^ 
cannot be explained by any of those standard Newtonian and Einsteinian dynamical effects 
which have been mismodeled (or unmodeled at all) in the force models of the EPM2008 
ephemerides. The same holds also for many exotic modifications of gravity proposed in 
the recent past to explain various kinds of cosmological/astrophysical observations. In 
particular, the DGP braneworld model is ruled out at about 3 — a level, while the existence 
of the force quadratic in the radial velocity proposed to explain the Pioneer anomaly must 
be excluded, at least in the spatial regions swept by the Saturn's orbit, at 11 — a level. 
Table [2] shows also that a finely-tuned cancelation of several unmodelled/mismodelled 
effects yielding just the estimated Arog^^t is unlikely. Both independent analyses by other 
teams of astronomers will be important in order to clarify the nature and the genuine 
existence of the anomalous behavior of the Saturnian perihelion. Moreover, it will be 
important for further studies to include the largest number of normal points as possible 
from spacecraft-based missions-in particular Cassini-covering the largest portion as possible 
of a full orbital revolution of Saturn. 
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Table 2: Nominal values, in " cy ^, of the secular precessions of wsat due to known classical and relativistic 
effects and by some non-standard forces; the effects with * have been modelled in EPM2008. The integrated 



N-body precession can be retrieved at http: / /ssd.jpl.nasa.gov/ txt /p-elem_tl .txt I have assumed for the 



quadrupole mass moment of the Sun J® = 2 x 10 ^ 



Sq = 190.0 X 10^9 kg m^ s-^ (IPiipers 



Fienga et al. 



Pireaux et al 



2007T ) and for its angular momentum 



19981 ). The precession by the the small asteroid ring has been computed 
(j2008f ) with TOiing = 0.34 x 10~^° Mq and n-ing = 2.8 AU. TNOs are the trans- 
Neptunian o bjects modelle d as a ring of mass totnOs and inner and outer radii i?inin and i?max, respectively 



2007H)- DGP is the 



(Dvali et al 



Lue and StarkmanI (|2003l ) perihelion precession in the Dvali-Gabadadzc-Porrati 
l2000f ) multidimensional brancworld model; tq w 5Gpc. Pioneer (Saturn) is the Pioneer anomaly 



at the Saturn's orbit ^Pio = —(1.8 ± 6.4) x 10 ^" m s (iNieto and Anderson 
Pioneer (w^) are the velocity-dependent forces proposed in (jStandish 



20051) . Pioneer {\vr\) and 
20081) : for them I used /C = 7.3 x 10"" 



s ^ and H = 6.07 x 10 m ^. The effect of a sp herically symmetric distribution of dark matter has been 



worked out in, e.g., (jKhriplovich and Pitieva 



20061 ) ■ while that of MOND is due to ( Sereno and Jetzer 



20061) 



with tmond = \/ GMq/Aq, Aq = 1.2 x 10 ^ ° m s ^. The apsidal precession induced by the cosmological 



constant A w 10 m ^ is due to 



Kerr et al. 



(2003|). The other parameters used are a, e,i, semimajor axis. 



eccentricity and inclination, respectively, of the planetary orbit, n — \J GMq /a^ Kcplcrian mean motion, 
Mq,Rq mass and equatorial radius, respectively, of the Sun, G Newtonian gravitational constant, c speed 
of light in vacuum. 



Dynamical effect 



Analytical expression 



w (" cy"i) 



N-body* 
Solar quadrupole* 
Small asteorid ring* 
TNOs* 
Schwarzschild* 
Lense-Thirring 

DGP 
Pioneer (Saturn) 
Pioneer {\vr\) 

Pioneer (w^) 
Dark matter 
MOND (fco = 1/2, m = 2) 
MOND (ko = 1, m = 1) 
Cosmological constant 



numerical integration 



3 _ 

2 ( 



■n.l2 ( Re \ ' 

l-e^)2 \—) 



1 



I sin^ i) 



Mq (-R„ax + flmi„)-R„ 

c^a(l-e^) 

4GSa 

c2o3(l-e2)3/2 

-p 3c 



. /a(l-£) 



2e-(l-e^)ln(l±f) 



7ina\/l — e^ ('—2 



47rGpdm\/l-e^ 
n 

(\ 2m 
rMOND J 

(\ 2m 
— 2— ) 
mOND / 



-1508.313 

3 X 10-^ 
1 X 10-5 

> 
0.014 

-1 X lO""^ 
TO.0005 
-12. 130± 43.130 
1.091 

-0.028 

> 
-1 X 10-5 

-8.098 

4 X 10-13 
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able 3: First row: esti mated correcti ons to th e standard pe rihelion pre cessions of Me rcury 



( iPitjeva 



2005r ). Venus flPitieva 



2008dh . Earth flPitieva 



2005|) and Mars flPitieva 



20051), in " 



cy ^. The quoted errors are realistic. Second row: anomalous perihelion precessions, in " 
cy~^, of Mercury, Venus, Earth and Mars induced by eq. (El). 



Mercury 



Venus 



Earth 



Mars 



-0.0036 ±0.0050 -0.0004 ±0.0005 -0.0002 ± 0.0004 0.0001 ± 0.0005 
-0.0011 -0.0017 -0.0019 -0.0023 
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